Ultra-high fidelity qubits for quantum computing 
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We analyze a system of fermionic 6 Li atoms in an optical trap, and show that an atom "on 
demand" can be prepared with ultra-high fidelity, exceeding 0.99998. This process can be scaled to 
many sites in parallel, providing a realistic method to initialize N qubits at ultra-high fidelity for 
quantum computing. We also show how efficient quantum gate operation can be implemented in 
this system, and how spatially resolved single-atom detection can be performed. 



The potential impact of quantum computing has stim- 
ulated a worldwide effort to developthe necessary exper- 
imental and theoretical resources [H, HJ . The most chal- 
lenging aspect of quantum computing is the requirement 
for ultra-high fidelity in each step. The necessary start- 
ing point is initialization of a scalable number of qubits 
at ultra-high fidelity. In this Letter, we consider a system 
of fermionic atoms and show that an atom "on demand" 
can be prepared in the ground state of a trap with fi- 
delity exceeding 0.99998, providing a realistic solution to 
the initialization problem in quantum computing. We 
discuss how that process can be scaled to many sites in 
parallel, enabling pairwise entanglement operations. Fi- 
nally we show how spatially resolved detection of each 
state can be implemented. 

We first address the question of how to trap one atom 
"on demand" in the ground state of a trap. Note that 
such a process should be scalable up to many optical 
traps in parallel. A single atom "on demand" in the 
ground state has not yet been demonstrated experimen- 
tally. However, progress has been made for bosonic atoms 
using the method of laser culling [3|, |j, [5[ . One of the key 
questions is that of fidelity of the number state, and in 
that regard bosons are not ideal because they rely on 
strong interactions to maintain a relatively large excita- 
tion gap, and to suppress low-frequency excitations dur- 
ing the culling process. This leads us to propose instead 
fermionic atoms where a precise number is rigorously en- 
forced by the Pauli exclusion principle. More specifically, 
we propose to use 6 Li as the atom of choice. This atom 
has the advantage that the interaction strength and sign 
(attractive or repulsive) can be tuned with an external 
magnetic field. Two magnetic sublevels of one hypcrfinc 
ground state, \F = \,mp = |) and \F = |, mp = — |) 
are used to define a qubit. We denote these states as ||) 
and ||) respectively. These states both become high- field 
seekers at large magnetic field with a well-defined fre- 
quency splitting that is nearly field-independent, hence 
insensitive to magnetic noise. 

The starting point of the "on demand" single atom 
preparation is laser-cooled 6 Li atoms that are optically 
trapped in the two spin states [a, uj. The atoms can 
be cooled by evaporation at a magnetic field around 300 



Gauss, where the scattering length a s w — 300ao (ao is 
the Bohr radius). a s is large enough for efficient evapo- 
ration of the spin mixture, and is also at a minimum as 
a function of magnetic field, greatly reducing the effects 
of magnetic noise. After evaporative cooling, a weakly 
interacting degenerate Fermi gas forms at temperature 
T<Tp, where Tp is the Fermi temperature. The single 
atom preparation process can be split into three steps: 

Step I The magnetic field is tuned to near ~ Gauss 
from the initial field of 300 Gauss, resulting in a 
non- interacting degenerate Fermi gas (DFG). In 
this state, a spin pair fills each level, up to approx- 
imately the Fermi level. 

Step II Atom pairs are ejected by laser culling. This is 
accomplished by adiabatic lowering of the optical 
potential. This prepares a single pair in the ground 
state. 

Step III The well is adiabatically split into two parts 
that are spatially separated. In the presence of a 
magnetic field bias this prepares one spin state on 
the left and the other on the right. Each atom can 
then serve as the initial state for a qubit. 

With this method, an array of N micro-traps would pre- 
pare 2N qubits. We now show in detail how ultra-high 
fidelity is enforced at each step. 

Step I In the presence of a scattering length a s w 
— 300a between the |f) and ||) states, the Fermi gas is 
weakly interacting. At T = 0, the Fermi gas may form a 
BCS state. The pairing gap for such a state can be es- 
timated through A O.bEp exp (-K/2kpa) Q. For a de- 
generate 6 Li gas with kp 1 ~ 1000oo, the pairing gap A w 
O.Q02Ep. With such a small A, the occupation probabil- 
ity of the lowest energy state (e.g. the k = state for a 
uniform gas) nk=o = (1 — ?7kAk) /2 ~ 1 — 4x 10~ 6 , where 

1 /2 

?7k = £k - E F , e k = h 2 k 2 /2m, e k = (A 2 + ry 2 .) is the 
quasi-particle excitation energy, and we take the chemical 
potential fi » Ep for the BCS state. In addition, for such 
a small pairing gap, finite temperature effects dominate, 
and there may even be no BCS pairing. However, the fi- 
nite temperature does not affect the ground state occupa- 
tion probability. Consider a temperature T = 0.05Tp, for 
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Figure 1: The simulation model for laser culling in a ID optical trap. We take xq = ^h/muj as the unit of length (m is the 
mass of the 6 Li atom), fko as the unit of energy, and hui/xo as the unit of force, a Truncated harmonic trap with magnetic 
gradient. The trapping potential for one state is shown. Parameters used to specify the trap: trap size, z, is the length of the 
parabolic portion in the potential profile; force, / is the slope of the linear portion in the potential profile. Dashed lines denote 
the levels of the ground and first-excited states, b Plot of P(E) as a measure of density of states of the quasi-bound states. 
Trap size z — 4.4, force / = 0.5. The vertical line (dotdashed) denotes the barrier height of the trap. Inset: Zoomed-in view 
of the resonance peak of the ground quasi-bound state (dots) as shown in the main figure and Lorentzian (solid) and Gaussian 
(dashed) regressions. Horizontal axis is in units of 7/2, where 7 is the FWHM of the resonance peak. 



which the ground state occupation probability is approx- 
imately 1/ [exp (-E F /k B T) + 1] = 1 - 4 x 1CT 5 . There- 
fore, a fast sweep of the magnetic field (i.e., the scattering 
length) to the non-interacting region does not affect the 
ground state occupation probability or the fidelity of the 
single atom preparation. 

Step II A non-interacting DFG in a deep optical trap 
serves as the starting point for the laser culling process. 
We assume that the trapping frequencies in the two trans- 
verse directions are significantly higher than along the 
axis, and only the ground transverse state can be sup- 
ported. We therefore limit our analysis to one spatial 
dimension. The trap wall is reduced to a level slightly 
below the ionization threshold of the first-excited state 
of the optical dipole trap. The trap reduction rate is cho- 
sen to fulfill the adiabaticity requirement. To that end, 
we maintain a constant trapping frequency, u>, through- 
out the laser culling process, which can be accomplished 
by dynamically varying the focus of the beam [9j. (In 
practice, this step may not be necessary, but it greatly 
simplifies the calculation.) According to the Wentzel- 
Kramcrs-Brillouin (WKB) method, atoms with energy 
much lower than the trap depth are completely unaf- 
fected by the change in the depth of the trap. The adi- 
abaticity condition is fulfilled as long as the WKB ap- 
proximation is valid, which holds until the trap depth is 
around 3huj/2, where h = h/2ir and h is Planck's con- 
stant. Near and beyond this point, the trap reduction 
rate must be slowed down to maintain adiabaticity. 

A constant force (or tilt), /, is also needed to sweep 
the atoms away from the micro-traps as soon as they are 
ionized, which is generated by applying a magnetic field 



gradient. For simplicity, we assume the force has a posi- 
tive sign for state ||). In our simulation, we approximate 
the optical trap with a truncated harmonic potential (see 
FigHh) which is specified by the trap size, z. After reach- 
ing a minimum trap size z rn , the trap is held for a certain 
time to allow ionized atoms to escape, while keeping the 
ground-state atom pair with high probability. The trap 
depth, magnetic force, and holding time are adjusted to 
optimize fidelity. Finally the trap wall (and size) is adia- 
batically raised to a higher level to preserve the resultant 
Fock state. 

Our simulation method is described below. We assume 
that the ground state fidelity is unchanged until the ion- 
ization energy of the first excited state is approached. 
Because of the tilt, the trap has no stationary bound 
state, only quasi-bound states, which, as z — > 00, be- 
come the bound states. The lifetimes of the quasi-bound 
states determine the rate of the change of trap occupa- 
tion probability. The optimized final trap depth, tilt, and 
holding time is calculated. 

Suppose a stream of incoming atoms is incident from 
x = —00, with energy E, scattered by the trap potential. 
Let tpE (x) be the wavefunction of the stationary state 
of the incoming atoms, which is assumed to have unit 
amplitude at x = — 00. We assume the trap is located at 
[—z/2,z/2]. Outside the trap (x < —z/2), tpE {%) is a su- 
perposition of Airy functions Ai(x) and Bi(x) with nec- 
essary superposition amplitudes and phases; inside the 
trap, 

4>E{x) = a(E)e-"'/ a H u (x), (1) 
where H v (x) is the Hermite function of degree v, a(E) 



is the amplitude. We require i/je (x) be continuous and 
differentiable for all x, by which a(E) is obtained. Of 
special interest to us are those states that have significant 
amplitude in the trap area, because they correspond to 
the quasi-bound states of the trap. We define P(E) = 
\a(E)\ . For simplicity, we can take P{E) as a measure 
of the density-of-states of the trap. We plot P(E) in Fig. 
[TJd. Note that in the limit of zero magnetic force, the 
peaks at E 0.366 and 1.29 correspond to the ground 
and first-excited states of the truncated harmonic trap, 
respectively. 

P(E) describes not only the ionization thresholds, but 
also the dynamical properties of the quasi-bound states. 
To see that, we study the evolution of a wavefunction 
(f>(x,t). Imagine that at time t = 0, we have 



<i>{x,0) 



0, 



-z/2 < x < z/2; 
otherwise, 



(2) 



where Eq is one of the resonance energies, c is a normal- 
ization factor such that \4>( x , 0) | 2 c?a; = 1 and there is 
no other atom source. Subsequently, the atom will start 
to tunnel out of the trap. The probability, R,E (t), for an 
atom to remain in the trap is given by f*^i 2 dx \(j>(x, t)| 2 . 
To evaluate Re (S), we expand <p(x,t) in terms of the 
wavefunctions of the stationary states, ipE^x). Note that 
in the vicinity of a resonance peak at E = Eq, one can 
write the wavefunction iPe(x) ~ C(E)a(E)ipE ( x ): where 
C(E) is a slowly varying quantity. Also 
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dxip E (x)4>E (x) 



dxip E ( x )^E (x) (3) 



for E ^ Eq, due to the oscillatory nature of Airy func- 
tions. We find 



R Eo (t) ^ \C(E )\' 



Eo+e/2 



P(E)e lEt dE 



Eo-e/2 



(4) 



where e is the range of integration. In the vicinity of a res- 
onance peak, the function P{E) is essentially Lorentzian 
(see inset of Fig. [TJd), and we finally obtain 



Re (t) 



-7En* 



(5) 



where -fE is the full width at half maximum (FWHM) 
of the resonance peak at E = Eq. The lifetime of the 
quasi-bound state at E — Eq is te — 7^ . Eq.© is 
used to determine an optimized combination of minimum 
trap depth, magnetic force, and holding time for the best 
fidelity. It is also worth noting that ultimately the dif- 
ference in the lifetimes between the ground and the first 
excited quasi-bound states determines the fidelity of pro- 
ducing a pair of atoms in the ground state of the trap. 

Taking into account all the steps in laser culling, we 
show the fidelity of preparing a single pair of atoms in the 
ground state of a micro-trap in Fig [2] Choosing culling 
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Figure 2: Ground state fidelity after laser culling. A varying 
holding time is allowed to ensure that the residual probabil- 
ity for atoms in excited quasi-bound states is no larger than 
10 -5 . Shown in the panel is the base-10 logarithm of the 
ground state fidelity loss. Red (blue) color represents high 
(low) fidelity. The white areas are out-of-range clippings: near 
the left-hand side, the white areas should be redder than its 
surrounding color; near center-top, the white area should be 
bluer than its surrounding color. 



parameters near the upper-left corner of the figure results 
in ultra-high fidelity. For a set of realistic parameters: 
trapping frequency to — 2tt x 1 kHz, magnetic field gra- 
dient 0.66 Gauss/cm, and truncated trap size 8.8 /mi, 
the ground-state to first-excited state lifetime ratio is 
7.53 x 10 5 . With a holding time of 218 ms, we get a 
residual probability of 10~ 5 for the excited-state and a 
ground state occupation probability larger than 0.99998. 

Step III Finally, the pair of atoms are split and 
trapped in the ground state of two adjacent micro-traps. 
In order to achieve deterministic splitting, where one spin 
state is driven to the left and the other to the right, 
we impose a magnetic field gradient (providing a force 
/) while we adiabatically split the optical tweezer (that 
holds the pair) into two beams [10|. With an appropri- 
ate magnetic field gradient and a separating barrier, this 
step can be realized at ultra- high fidelity. Since |J.) is a 
low- field-seeker and |f) is a high-field-seeker at low mag- 
netic field, each atom is displaced to a different location 
as soon as the trap is split. 

We adopt a simpler potential than that of a realistic 
optical tweezer. The trap is composed of two spliced 
sections of parabolic trap with identical trapping fre- 
quency lu, each of which has an energy minimum, lo- 
cated at ±d/2, respectively, where d is the separation 
distance. We numerically calculate the eigen-energy lev- 
els and wavefunctions of such a double- well structure in 
the parameter space of / and d. One needs to choose 
a path in the parameter space such that a sufficient gap 
between the ground and first-excited states is maintained 
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throughout the splitting process in order to suppress 
transitions from the ground state. With a trapping fre- 
quency uj = 2n x 1 kHz, a magnetic field gradient of 0.66 
Gauss/cm and a separation displacement d = 6.25 /mi, 
we find 0.99998 splitting fidelity. 

So far, we have shown that two initialized fiducial 
states can be prepared at ultra-high fidelity. How can 
one scale this process to multiple sites, and make the 
"switchyard" of multiplexed beams to perform the re- 
quired complex operations? The scalability of this sys- 
tem can be achieved using scalable microelectromechan- 
ical systems (MEMS) technology Using this tech- 

nology, an array of beams can be directed to each site 
where individual atoms can be trapped. An alternative 
approach was developed in ref. [l2|. By steering the 
beams, one can transport individual atoms and bring 
them into pairwise interaction with arbitrary control. An 
array of micro-traps would most easily be accomplished 
with red-detuned beams that create attractive potentials 
along the axis. This technique also enables the creation 
of a two-dimensional optical trap array, and entangle- 
ment of any pair using the qubit transfer technique of 
ref. [l(J. This could extend the linear case where only 
nearest-neighbor operations are possible. 

We now consider single-qubit gates and two-qubit 
gates in this system. The implementation of a single- 
qubit gate requires the capability to address each atom 
individually. This can be accomplished with stimulated 
Raman transitions, as is currently employed with trapped 
ions [HI, [3] and with neutral atoms (lH . The realization 
of a two-qubit gate can be based on collisions between 
bosonic atoms [16( ■ This scheme was realized experimen- 
tally with atoms in an optical lattice 17). Fermionic 



atoms in the same internal state cannot collide due to 
the Pauli exclusion principle. However, atoms in differ- 
ent states can have a large collisional shift, which can be 
used to engineer two-qubit SWAP (or v / SWAP) gates, as 
proposed in ref. 18j . The scattering length can be made 
very large by tuning closer to the Feshbach resonance. 
For a set of parameters uj = 2?t x 1 kHz and a s ~ 330ao, 
a high fidelity VSWAP gate can be implemented in a time 
period of about 40 ms. We envision that with an array of 
many qubits, a sequence of VSWAP operations can build 
scalable entanglement in this system. This two-qubit op- 
eration, together with single-qubit rotations, provides a 
set of universal quantum gates. 

The detection of each qubit state at the end of a 
quantum computation can be accomplished by spatially- 
resolved fluorescent imaging This stage must 
also resolve the spin of each location. This can be accom- 
plished with the same method that was used to separate 
the spin pair, followed by fluorescent imaging. 

One important consideration is the limitation on fi- 
delity set by decoherence. One fundamental source is 
light scattering by the optical tweezers. The scattering 
rate is proportional to the third power of the frequency. 



This was the motivation for using a CO2 laser for opti- 
cal trapping 0]. In our case, such a laser can be used 
for trapping in two transverse directions where the trap 
depth must be higher. This would confine the atoms 
along a line. However, the optical tweezers that imple- 
ment logic gates must be at a shorter wavelength in order 
to have sufficient spatial resolution. The trap depth in 
this direction need only be enough to contain the ground 
state, so we estimate that light scattering will not limit 
the quoted fidelity. Other sources of decoherence are laser 
intensity noise and pointing noise. The fundamental limit 
in this case is shot noise, and should also not limit the 
predicted fidelity. 

In conclusion, we find that fermionic atoms can be used 
to realize ultra-high fidelity for quantum computing. Fu- 
ture work will analyze the optimum strategies for pair- 
wise entanglement in this system. 
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